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摘要 : 对 大 型 射电 望远镜 主 反射 面 的 单 块 面板 变形 进行 了 研究 ， 将 机 器 人 学 六 自由 度 姿 态 变 化 
应 用 至 大 型 射电 望远镜 反射 面 面 板 的 变形 研究 当中 ， 提 出 了 一 种 基于 刚体 六 自由 度 变形 的 面板 
变形 表征 方式 。 基 于 该 刚体 六 自由 度 姿态 变形 理论 ， 提 出 了 一 种 基于 光学 方法 的 天 线 面板 变形 重 
构 技 术 ， 该 项 技术 通过 固 结 于 面板 上 的 激光 单元 所 呈现 的 光斑 图 像 来 快速 高 精度 地 重 构 出 面板 
的 变形 。 深 入 地 探究 了 激光 光斑 与 天 线 面板 变形 两 者 之 间 的 映射 关系 ， 将 天 线 面板 空间 变形 问题 
转化 为 有 限 区 域内 的 激光 光斑 位 移 问题 ， 对 该 方案 的 光斑 图 像 进行 了 研究 ， 提 供 了 有 效 的 光斑 医 
像 处 理 算法 ， 将 图 像 差 分 方法 与 Canny 边缘 检测 算法 相 结合 ， 能 够 有 效 地 解决 传统 Canny 算法 
在 处 理 小 目标 时 的 问题 ， 还 采用 了 高 精度 亚 像素 质心 定位 算法 ， 能 够 对 光斑 点 进行 亚 像素 级 别 的 
定位 。 为 天 线 面板 的 测量 提供 了 一 种 新 的 、 可 靠 的 图 像 处 理 方法 ， 也 为 提高 天 线 面板 测量 的 精度 
和 效率 提供 了 新 思路 和 方法 。 
关 键 词 : 大 型 射电 望远镜 ;变形 重 构 ， 图 像 处 理 ， 亚 像素 
中 图 分 类 号 : P111.44 文献 标识 码 : A 


1 引 言 


随 着 各 项 深 空 探测 活动 如 宇宙 上 暗物质 探索 、 宇 宙 演 化 研究 、 射 电 天 体 探索 等 的 开展 ， 天 
线 作为 接收 和 传输 信号 的 重要 工具 ， 也 不 断 朝 着 口径 更 大 和 频率 更 高 的 方向 发 展 。 天 线 的 
径 越 大 ， 在 观测 深 空 物体 时 能 够 获得 的 信 噪 比 越 高 ， 及 得 到 的 观测 结果 越 准确 。 但 是 ， 天 线 
口径 越 大 其 表面 精度 就 越 难以 维持 ， 这 意味 着 对 天 线 表面 精度 的 要 求 也 随 着 口径 的 增 大 而 
变 得 更 加 苛刻 。 同 时 ， 随 着 频率 的 提高 ， 对 于 天 线 的 反射 面 面 型 精度 也 有 更 高 的 要 求 。 在 高 
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频率 下 ， 即 使 微小 的 表面 形变 也 会 对 反射 产 4 
需要 使 用 更 高 精度 的 测量 系统 ， 对 天 线 的 面 型 


精度 。 


天 文学 进展 


42 卷 


E 明 显 的 影响 ， 导 致 信号 的 损失 或 扭曲 。 因 此 ， 
进行 及 时 校正 和 补偿 ， 以 确保 天 线 表面 的 形态 


EAE € DESA : 


Ruze 公式 给 出 了 天 线 增益 G 与 其 主 反射 面 表面 

G 4TE 

gp - x 

Heb, Go 是 理想 的 、 无 误差 的 反射 面 增益 ，G 表示 天 线 增益 ， 和 表示 波长 ， 单 位 为 米 (m); 


= 表示 天 线 表面 精度 ， 通 常 以 米 (m) 或 毫米 (mm) 为 单位 。 此 式 也 再 次 说 明了 要 想 保证 天 线 
工作 效率 必须 长 期 保持 较 高 的 天 线 表面 精度 。 


) | (1) 


现 有 的 天 线 测量 系统 在 一 定 程度 


大 型 天 线 或 高 频段 上 需要 更 高 精度 的 


上 受 到 测量 速度 、 测 量 精度 以 及 设备 的 限制 ， 特 别 是 在 
测量 方法 和 设备 。 然 而 ， 有 些 测量 系统 需要 进行 网 格 
划分 或 点 扫描 等 较 耗 时 的 操作 ， 因 此 测量 速度 较 慢 。 如 相位 全 息 法 ”需要 进行 远 场 的 幅 值 和 


相位 测量 ， 对 于 大 型 天 线 来 说 需要 测量 大 量 的 数据 点 ， 导 致 测量 时 间 较 长 ， 并 且 全 居 法 受信 
号 源 限制 ， 很 难 实现 全 方位 、 任 意 俯仰 角 下 的 测量 。 


~ 


要 求 ， 只 有 在 较 好 的 天 气 条 件 才 能 获得 准确 数据 ， 测 量 精度 很 难 达到 0.2 mm 以 内 。 此 外 ， 
有 些 测 量 系 统 还 需要 一 些 昂贵 的 测量 设备 ， 如 激光 雷 
高 ， 这 就 限制 了 测量 设备 的 使 用 范围 和 普及 度 。 摄 
量 精度 受 天 线 口 径 的 影响 很 大 ， 口 径 越 大 测量 精度 越 差 ， 因 此 仅 用 在 面 形 误差 要 求 较 低 或 


E 


又 如 激光 测 距 法 ”还 对 环境 条 件 有 一 定 


达 扫 描 仪 和 激光 跟踪 仪 所 等 ， 成 本 非常 


HU, 
H2 


测量 法 的 测量 过 程 相对 简单 ， 但 其 测 


径 较 小 的 天 线 上 。 综 上 所 述 ， 现 有 测量 


方法 在 大 型 射电 望远镜 实时 测量 方面 仍然 存在 测量 时 


间 较 长 、 测 量 过 程 繁 琐 、 易 受 环境 影 


此 对 面板 进行 研究 是 必要 的 。 本 文 基 于 大 型 双 反 射 开 


犁 、 成 本 高 等 不 足 之 处 。 
射电 望远镜 主 反射 面 的 面 型 变化 是 一 个 复杂 的 问题 ， 由 于 反射 面 是 由 多 个 面板 组 成 ， 因 


射电 望远镜 天马 望远镜 号 的 主 反射 


完成 单 块 面板 的 测量 。 这 意味 着 无 论 


将 天 线 面板 空间 变形 问题 转化 为 有 限 


参考 。 


面 面 板 展开 了 变形 研究 ， 提 出 了 一 种 可 以 实时 测量 面板 变形 的 光学 映射 测量 方法 。 该 方法 无 
需 任何 昂贵 的 设备 以 及 耗 时 的 操作 ， 只 需 利用 相机 i 


己 录 光斑 位 置 ， 然 后 由 论文 算法 即 可 快速 


_ 


企 系 统 复 杂 度 上 还 是 在 测量 成 本 和 测量 时 间 上 都 将 被 
大 大 缩小 ， 此 外 该 测量 方法 不 受 天 线 姿态 和 仰角 以 及 馈 源 的 限制 。 本 文 提 出 的 光学 映射 方法 
区 域内 激光 光斑 移动 问题 ， 这 一 创新 解决 了 大 口径 射 
电 望 远 镜面 板 变形 测量 中 的 诸多 难题 ， 并 为 不 断 探索 和 创新 面 形 测量 新 技术 提供 了 有 益 的 


2 基于 光学 映射 方法 的 面板 变形 重 构 技 术 


本 文 的 面板 变形 重 构 技 术 基 于 光学 映射 方法 ， 包 括 三 个 核心 部 分 : 刚体 六 自由 度 理论 、 


面板 变形 与 光斑 映射 关系 的 数学 模型 


光斑 位 置 的 动态 变化 来 实时 映射 面板 


以 及 图 像 处 理 技术 。 该 方法 通过 空间 六 自由 度 描 述 和 
记录 刚体 面板 的 位 置 和 姿态 变化 ， 即 刚体 变形 ， 然 后 建立 面板 与 光斑 之 间 的 映射 关系 ， 利 用 
在 空间 中 的 6 个 自由 度 变化 。 相 机 作为 视觉 系统 拍摄 


1 期 BRE, T: 基于 激光 映射 方法 的 天 线 画 


光斑 位 置 ， 并 通过 本 文 的 图 像 处 理 技术 提取 光斑 信 


的 空间 六 自由 度 变 量 。 


本 文 主要 采用 光学 映射 方法 的 天 线 面 


j 板 变形 重 构 技 术 研究 
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最 终 利用 光斑 位 置信 息 计 算 刚体 面板 


变形 测量 技术 ， 原 理 如 图 1 所 示 。 首 先 在 天 线 面 


板 4 个 角 点 上 分 别 设置 激光 器 ， 并 利用 激光 接收 面板 接收 激光 器 发 出 的 射线 " ， 形 成 二 维 的 
光斑 图 像 ， 并 通过 一 台 高 像素 工业 相机 作为 拍摄 系统 "记录 光斑 的 位 移 。 当 天 线 面 板 未 发 生 
变形 时 ， 每 个 射线 与 接收 面板 的 交点 是 均匀 分 布 的 ， 当 天 线 面 板 发 生变 形 时 ， 每 个 射线 与 接 


收 面板 交点 的 光斑 位 置 也 会 发 生 改变 。 
板 上 点 的 变化 重 构 面 板 变 形 。 基 于 上 述 思 


通过 推导 出 两 者 之 间 的 映射 模型 ， 
合 射电 望远镜 的 整体 结构 ， 本 文 以 天 马 望 远 


便 可 以 依据 接收 面 


镜 主 反 射 面 的 单 块 面板 作为 研究 对 象 ， 开 展 了 光学 映射 方法 的 天 线 面板 变形 重 构 技 术 研究 。 


图 1 光学 映射 方法 天 线 面板 变形 测量 原理 图 


3 面板 变形 与 光斑 映射 关系 的 数学 模型 


3.1 “六 自由 度 参数 表征 及 坐标 系 建立 


方向 上 的 轴 进 行 旋转 ， 因 此 一 共有 6 个 独立 的 自由 度 
可 以 实现 在 三 维 空间 中 的 各 种 运动 ， 例 如 在 三 维 


里 ， 假 设 天 线 面 板 为 刚体 ， 便 可 以 使 用 


的 六 自由 度 变 换 和 


间 中 的 位 置 和 姿态 ， 就 可 以 表征 天 线 面 板 的 变形 。 


机 器 人 学 中 的 “自由 度 ”(degree of freedom, DOF) 指 的 是 机 器 人 在 运动 过 程 中 独立 的 
运动 自由 度 。 一 个 机 器 人 的 DOF 数量 决定 了 其 在 空 
由 度 运动 是 指 机 器 人 可 以 在 三 维 空间 中 治 着 zx、 从 z 这 


zr E 够 执行 的 运动 种 类 和 数量 。 六 自 


这 3 个 方向 平移 ， 同 时 可 以 绕 这 3 个 
。 通 过 这 6 个 自由 度 的 组 合 ， 机 器 人 


+ 抓 取 物 体 、 对 物体 进行 加 工 等 。 同 
E 阵 来 描述 刚体 在 三 维 空 


N 


对 天 线 面 板 不 同 姿态 进行 表征 时 ，6 个 自由 度 参数 分 别 是 沿 着 六 y. z 这 3 个 轴 向 的 位 


FER Ar. Ay. Az, DUREE. y. “这 3 个 轴 向 的 人 
手段 我 们 可 以 反 疝 习 


征 主 面 的 面板 变形 ， 并 且 利 用 这 一 
其 在 不 同 姿态 下 的 形变 情况 。 


为 了 更 好 地 描述 天 线 面板 的 变形 ， 我 们 针对 单 


Vou 


E a、B、y。 利 用 六 自由 度 参 数 表 
E 构 天 线 面 板 的 变形 ， 并 准确 地 测量 出 


四 板 建 立 了 坐标 系 ， 便 于 对 其 进行 分 
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fT 和 建 模 。 具 体 如 图 2 所 示 ， 以 面板 的 中 心 点 为 坐标 原点 o， 过 中 点 且 垂直 于 面板 的 轴线 为 
E 标 轴 z， 以 面板 平面 为 z-o-y 面 建 立 单 块 面板 的 单元 坐标 系 。 


DE eb 


| 单元 坐标 系 
EU CT 
2 | 


全 局 坐标 系 


2 ”方案 设计 与 面板 单元 坐标 系 


32 ” 主 面 变形 - 副 面 光斑 方程 的 建立 

对 于 单 面板 来 说 ， 仅 需 构 建 该 面板 的 位 移 与 该 面板 上 激光 发 射 单元 投影 在 副 面 上 
的 激光 光斑 位 置 变化 之 间 的 映射 关系 即 可 。 对 于 刚体 面板 ， 考 虑 其 空间 六 自由 度 位 移 
(Ac, Ay, Az, o, B, 7)， 根 据 刚体 平 移 旋 转变 换 和 矩阵 ， 可 以 将 其 变形 前 后 的 坐标 变换 关系 
表示 为 : 


Ta To 
“= T(Az, Ay, Az, a, B, y) x Yo ， (2) 
Zd Zo 
1 1 


其 中 ，(zo, yo, zo) 表示 面板 变形 前 刚体 上 固 结 点 在 单元 坐标 系 下 的 坐标 ，(ze， ya, za) 表示 
面板 变形 后 刚体 上 固 结 点 在 单元 坐标 系 下 的 坐标 。 


T(Az, Ay, Az, Q, B, y) = T4 X Ts x 了 2 x Ti , (3) 
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0 0 0 
cosa 一 Sna 0 
sina cosa 0 

0 0 1 


[em Mem e E 


cosa 0 sing 


0 
0 1 0 0 
T, = 
一 Sn 0 cos8 0 
1 


0 0 0 
, (4) 
cosy —siny 0 0 
T, — siny cosy 0 0 
0 0 1 0 
0 0 0 1 
100 Ar 
010 A 
Ti 4 
0 0 1 Az 
000 1 


Hop. T(Az, Ay, Az, a, B, vy) Sez WA ETE EREE, Ti. To. Ts 表示 三 个 轴 向 的 
旋转 和 矩阵，Ta 表示 平移 矩阵 。 
3.3 ” 非 线 性 方程 组 最 优化 求解 问题 

将 面板 4 个 角 点 发 出 的 激光 束 视 为 一 个 刚性 杆 ( 见 图 2)， 因 其 与 刚性 面板 固 结 ， 所 以 
面板 和 激光 束 可 视 为 一 个 整体 ， 即 跟随 面板 发 生 六 自由 度 位 移 。 所 以 无 论 是 面板 上 还 是 激 
光束 上 ， 变 形 前 后 的 点 均 可 以 用 式 (2) 来 表示 。 空 间 中 线段 与 曲面 相交 ， 则 曲面 上 的 交点 
(x, y, z) 应 满足 : 


TX VY XT (5) 
T2 — Tı Y2 — yı Z2 — 21 


其 中 ，(zi, y, 21) (z2, Y2, 22) 分 别 表示 线段 上 的 两 点 ， 在 本 文中 为 线段 起 点 (面板 的 角 
点 ) 以 及 线段 终点 。 
天 线 副 反射 面 在 全 局 坐标 系 下 的 曲面 方程 为 : 


11 
za y ry)rlcl-a ， (6) 


RF, as 表示 副 反 射 面 的 实 半 轴 ，bs 表示 副 反 射 面 的 虚 半 轴 ，! 表示 主 反射 面 顶点 与 副 反 射 
面 顶 点 之 间 的 距离 ，z、y、z 表示 副 反射 面 上 的 点 在 全 局 坐标 系 中 的 空间 坐标 。 
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XX (5) 可 视 为 一 个 方程 组 ， 对 于 单个 面板 ， 可 列 出 4 个 方程 组 : 


HF, (£m, Ym, zm) AED 
F 的 坐标 (本 文 取 F 为 光束 与 副 面 的 交点 ， 与 M 点 不 同 的 是 变形 后 
量 的 参数 形式 )，(zi, yi, Zi) 表示 面板 上 4 个 角 点 的 坐标 。 


Xm — Xf Z Um — Yf 


Ui — Yf 
Xm — Tf - £m — £f 


(i=1, 2,3,4) , 


后 光斑 点 M 的 中 心 4 


Zi — Ef 
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(T) 


EER, (vp, yp, zt) 表示 激光 束 上 的 一 点 
HP 坐 标 是 含有 待 求 变 


面板 发 生变 形 前 所 有 坐标 是 已 知 的 ;发 生变 形 后 ， 根 据 亚 像素 算法 获取 光斑 质心 坐标 


他 坐标 为 待 求 坐 标 。 


M (zm, Ym, Zm)» 2 


示 为 已 知 的 变形 前 多 
参数 (Ax, Ay, Az, 


计算 得 到 这 6 个 参数 ， 便 很 容易 对 本 
角 点 都 可 以 列 出 一 个 形 如 式 (7) 的 方程 组 ， 包 含 两 个 方程 。 那 么 对 于 6 个 变量 ， 


Q, B, y) 


民 据 式 (2)， 可 将 


标 与 平移 旋转 变换 矩 


意 3 个 角 点 便 可 实现 静 定 问题 (6 个 方程 、 
点 得 到 的 静 定 解 并 不 一 定 满足 所 有 方程 ， 即 


等 式 均 被 满足 ， 可 认为 存在 误差 的 情况 下 ， 方 程 组 为 不 存在 解 的 矛盾 方程 组 。 
化 为 最 优化 问题 ， 即 需要 找到 一 组 六 自 
相差 最 小 。 本 文采 用 梯度 下 降 算 法 


解 ， 得 到 式 (7) 的 最 优 解 。 


4 图 像 处 理 算 法 


由 上 述 可 知 ， 天 线 面 板 的 变形 是 通过 光斑 位 置 变化 来 重 构 的 ， 因 
究 过 程 中 我 们 发 现 ， 光 斑 的 最 小 位 移 量 在 微 
立方 法 完成 。 


移 才 能 重 构 精 确 的 面板 变 
立 移 


米 级 ， 而 微米 级 的 


由 度 


E. (UK Bn 


(£f, Yf, zp) 以 及 (£i, Yi, zi) * 
阵 相 乘 ， 可 将 待 求 参 数 转 变 为 六 自由 度 平移 旋转 


板 的 变形 进行 数学 表征 。 对 于 单 块 面板 ， 每 一 个 


只 需 其 中 任 

6 个 未 知 数 ) 的 求解 。 但 考虑 实际 情况 下 3 个 角 

组 六 自由 度 平 移 旋转 参数 无 法 使 方程 组 每 个 

因此 将 问题 转 

F 移 旋转 参数 ， 使 得 式 (7) 中 所 有 方程 等 式 两 边 

: (Levenberg-Marquardt) 对 上 述 非 线 性 最 优化 问题 进行 求 
此 需要 精确 的 光斑 位 


但 事实 上 视觉 系统 托 


摄 的 图 


Canny 边缘 检测 等 


4.4 图像 差 分 法 


Z 


以 用 于 检测 两 幅 图 


图 像 


RZ ERRE. H 
值 ， 从 而 检测 目标 是 否 发 生 了 变化 。 
或 者 出 现 或 者 消失 等 引起 的 。 因 
稳定 的 情况 下 ， 差 分 法 可 以 很 好 地 完成 这 个 任务 。 针 对 本 方案 ， 


像 差 分 法 (image difference method)" 是 一 


使 得 视觉 系统 无 法 轻易 识别 ， 


因此 需要 借助 亚 像 素 定 


像 是 存在 噪声 以 及 复杂 
预 处 理 方法 。 这 种 方法 在 处 理 较 小 的 光斑 时 十 分 有 效 ， 因 


因此 本 文 引 入 了 


背景 的 ， 


,基于 图 


基本 原理 是 将 两 幅 图 像 ; 


It, B 


通常 1 


KHMER H 


图 像 差分 法 和 


为 光斑 的 
像素 尺寸 较 小 ， 对 于 微小 的 噪音 更 为 敏感 ， 而 边缘 检测 可 以 更 精确 地 提取 光斑 的 轮廓 信息 ， 
从 而 减少 噪音 对 质心 定位 的 影响 ， 实 现 光斑 质心 坐标 的 实时 精准 提取 ”。 


的 ， 需 要 先 拍 摄 


像 灰 度 差异 的 目标 检测 方法 ， 可 
行 相 减 ， 得 到 像素 级 别 的 差异 


况 下 ， 两 幅 图 像 之 间 的 差异 是 由 目标 物体 的 移动 
像 差分 法 常 被 用 于 消除 目标 图 像 的 背景 。 在 背景 相对 
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无 目标 物体 的 纯 背 景 图 作为 参考 图 像 ， 然 后 将 含 目标 物体 的 图 像 与 参考 图 像 进行 差分 ， 得 到 
目标 物体 的 位 置信 息 图 像 ， 从 而 实现 背景 消除 的 效果 。 

差分 方法 应 用 于 本 方案 中 ， 可 以 很 好 地 消除 图 片上 不 相关 的 复杂 背景 以 及 强 光 源 的 干 
扰 。 对 于 本 方案 的 其 体 实现 过 程 如 图 3 所 示 : 在 激光 器 未 开启 之 前 先 拍 摄 无 激光 光斑 的 原 
始 背景 图 像 b)， 然 后 打开 激光 器 拍摄 光斑 图 像 a)， 最 终 用 a) 减 去 b) 便 可 得 到 去 除 背 景 医 
像 的 激光 光斑 图 像 c)。 可 见 差分 法 对 背景 以 及 强 光 源 的 消除 效果 是 非常 理想 的 。 


H 标 激 治 点 不 相 zs dj m 


图 3 图像 差分 法 


4.2 ”光斑 图 像 边缘 处 理 算法 

上 述 图 像 差分 法 可 以 较 好 地 滤 除 复杂 的 背景 以 及 一 些 强 光源 的 干扰 ， 但 是 不 可 避免 地 
一 些 弱 光 依然 存在 。 这 一 情况 主要 是 由 参考 图 像 与 目标 图 像 拍 摄 的 时 间 差 导致 ， 即 无 法 拍 
摄 出 与 目标 图 像 完全 相同 的 背景 参考 图 。 这 些微 小 噪声 一 般 比 目标 光斑 的 灰 度 值 低 (亮度 更 
低 )， 有 的 甚至 肉眼 无 法 识别 ， 但 是 依旧 会 对 后 续 亚 像素 级 别 的 质心 定位 造成 影响 。 边 缘 检 
测 技术 可 以 精准 有 效 识别 光斑 的 边缘 轮廓 ， 并 达到 滤 除 弱 光 、 增 强 光斑 的 效果 。Canny 边 
缘 检测 “是 目前 应 用 最 为 广泛 的 一 类 边缘 算法 。 传 统 的 Canny 边缘 检测 流程 如 图 4 所 示 ， 
主要 分 为 5 部 分 : 高 斯 滤波 、 像 素 梯度 计算 、 非 极 大 值 抑制 、 双 阔 值 检测 和 连接 边缘 。 


图 像 输 入 高 斯 滤波 像素 梯度 与 方向 计算 FFRAE SRRA ”边缘 连接 


图 4 Canny 边缘 检测 流程 示意 图 


传统 的 Canny 算法 在 处 理 较 大 目标 主体 时 有 很 好 的 识别 效果 ， 但 是 对 于 本 方案 中 的 激 
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光 光 斑点 这 样 的 小 
问题 ， 本 文 提出 一 种 新 的 光斑 


目标 ， 传 统 


的 Canny 算法 容易 受到 
图 像 处 理 思路 ， 将 


天 文学 进展 


Ab Ed Ri 


啊 而 难 


H m 


图 


合 。 在 图 像 预 处 理 阶段 ， 采 
点 更 加 突出 ， 然 后 再 使 用 Canny 
db =E 


MEG 


p 


f 


YA 


1 


j 图 像 差分 方法 对 原始 
了 边缘 检测 。 


Z] 


H 


a 


ER 


PRA 
像 差 分 方法 与 Canny 边缘 检测 算 ; 
像 进行 处 理 ， 去 除 背 景 干扰 ， 


| 


LI 


通过 这 种 


Ti 
F 


斑点 准确 识别 的 同时 ， 有 效 去 除 背 景 
的 实验 结果 可 以 看 出 ， 传 统 的 Canny 算 


扰 ， 提 高 光斑 质心 的 准 在 
法 虽然 在 大 目标 上 有 较 理 


中 合 的 算法 ， 


Hj 


性 


原始 图 像 


接应 用 于 本 方案 ， 本 文 提出 的 融合 算法 更 加 适 


用 于 本 方案 的 激 》 


6 光斑 图 


和 a 
相 
JO 


的 处 理 


"T 
定性 。 对 比 


结果 ， 


像 的 边缘 检测 。 


传统 Canny 处 理 


但 是 


原始 图 像 AK BE RE RE 密度 0.02 椒 盐 噪声 边缘 处 理 结 果 
5 ”基于 小 目标 检测 的 Canny 边缘 算法 
4.3 ” 亚 像素 光斑 质心 定位 
图 像 经 上 述 处 理 之 后 便 已 完成 图 像 的 预 处 理工 作 ， 接 下 来 只 需 将 处 理 后 的 仅 含 有 光斑 位 
置信 息 的 图 像 进 行 质心 定位 算法 便 可 完成 对 光斑 的 质心 提取 。 同 样 ， 亚 像素 质心 算法 “也 
是 一 种 基于 图 像 灰 度 值 的 计算 方法 ， 但 是 区 别 于 差分 法 和 边缘 检测 等 整 像素 级 别 的 算法 ， 亚 
像素 质心 算法 则 是 基于 


S 


灰 度 中 心 法 是 最 


M «AN 
u 2264 22i Tij 


Te 


像素 级 别 ， 有 更 高 的 计 
简单 的 种 亚 像 素质 心算 法 ， 计算 过 程 为 : 


精度 。 


M CN 
i 25-4 2d Yij 


C 


MxN i 
其 中 ，ze 和 y 分 别 表示 亚 像 素质 心 坐标 ;AM 为 
yi 为 光斑 中 点 的 像素 坐标 。 利 用 


MxN ' 


图 像 像素 高 度 ， 为 


(8) 


图 像 像素 宽度 ， zx;; 和 
该 公式 遍历 图 像 的 每 一 个 像素 点 ， 计 算 其 位 置 坐标 以 及 像 
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素 值 得 到 两 个 变量 zi 和 y;;}， 并 将 其 坐标 与 像素 值 累 加 。 遍 历 完整 个 图 像 后 ， 用 累加 后 的 
结果 分 别 除 以 图 像 所 有 像素 ， 得 到 的 结果 即 为 光斑 质心 的 横 坐 标 和 纵 坐 标 。 
重心 法 也 可 以 视 为 灰 度 质心 法 的 一 种 ， 重 心 法 对 于 图 像 的 要 求 一 般 比 较 低 ， 适 用 于 多 种 
类 型 的 图 像 处 理 任务 。 在 重心 法 中 ， 像 素 的 灰 度 值 作 为 权重 ， 用 于 计算 光斑 的 质心 坐标 : 


PS mu pe Tij X G(x, y)? uus na NS yij X G(z,y)^ 
x m G(z, y)* l pue um G(v,y)* 


其 中 ，ze 和 w 表示 亚 像素 质心 坐标 ，G(z,， y) 为 像素 点 的 灰 度 值 ，a 为 权重 系数 ， 权 重 系 
数 越 大 ， 则 亮度 高 的 像素 点 权重 越 大 ，z;; 和 y;; 为 光斑 中 点 的 像素 坐标 ;JM 为 图 像 像素 高 
E, N 为 图 像 像素 宽度 。 在 计算 质心 坐标 时 ， 灰 度 值 更 高 的 像素 点 会 被 赋予 更 大 的 权重 ， 
因此 亮度 更 高 的 像素 点 对 于 质心 坐标 的 计算 会 有 更 大 的 贡献 。 

海 森 法 是 一 种 寻找 函数 极 值 的 优化 算法 ， 可 以 利用 海 森 矩阵 确定 光斑 灰 度 函数 极 值 ， 从 
而 得 到 光斑 中 心 的 亚 像素 位 置 。 在 这 种 方法 中 ， 海 森 矩 阵 用 来 计算 目标 函数 的 二 阶 偏 导数 ， 
以 确定 目标 函数 在 该 点 的 极 值 位 置 。 根 据 海 森 法 的 原理 ， 目 标 函 数 极 值 点 的 方向 是 海 森 矩 阵 
的 特征 向 量 ， 而 极 值 点 的 值 则 是 海 森 矩阵 特征 值 的 相反 数 。 因 此 ， 在 确定 光斑 中 心 的 亚 像 素 
位 置 时 ， 通 常 使 用 海 森 矩 阵 的 特征 向 量 方向 来 指导 质心 计算 的 方向 ， 进 而 得 到 激光 光斑 中 心 
的 亚 像素 位 置 。 知 定义 图 像 灰 度 函 数 为 f(z1， zx2， Zz3,… ，zn)， 则 可 以 定义 函数 了 的 海 森 
AREE: 


(9) 


9? f O0? f 
O271 i OZ17Zn 
H(f) = : En : i (10) 
of of 
Ln T1 v O27Tn 
高 斯 拟 合法 求 质心 坐标 是 一 种 基于 高 斯 拟 合 的 亚 像素 定位 算法 。 该 方法 的 基本 思想 是 


将 光斑 看 作 是 二 维 高 斯 分 布 ， 将 像素 值 看 作 高 斯 分 布 函数 的 采样 值 ， 利 用 最 小 二 乘法 拟 合 高 
斯 分 布 函数 ， 得 到 拟 合 参数 ， 包 括 高 斯 分 布 的 中 心 、 标 准 差 和 幅 值 等 。 根 据 拟 合 参数 计算 光 
斑 的 质心 坐标 ， 其 中 光斑 的 中 心 即 为 高 斯 分 布 的 中 心 ， 具 体 实 现 过 程 为 : 


Ei ligje (i70) /20° n l;y;e- (iv /20* 


Uo l;e- (i7 20)? /20? ， Yc 二 X l;e- Wi- v9)? /20? , (11) 


Te = 


其 中 ，ze 和 yy 是 光斑 的 质心 坐标 ，z; 和 y; 是 第 i 个 像素 的 坐标 ，1; 是 第 i 个 像素 的 灰 度 
值 ，zo 和 wo 是 光斑 中 心 的 估计 坐标 ，o 表示 二 维 高 斯 分 布 函数 的 标准 差 ，n 表示 采样 点 的 
数量 ， 也 即 图 像 中 用 于 计算 光斑 质心 坐标 的 像素 数量 。 

4.4 ”光斑 质心 坐标 求解 

通过 上 述 图 像 处 理 以 及 质心 定位 算法 ， 便 可 以 得 到 目标 光斑 在 图 像 上 的 二 维 质心 位 置 
坐标 C(xe, yc) 


Ces Ye) L7 


f 
Sg) D’ MS 
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其 中 ，C(ze yc) 表示 光斑 中 心 在 相机 像素 平面 中 的 坐标 ;，S(x, y) 表示 副 反 射 面 上 的 激光 光 


斑 的 二 维 (x, y) 坐标 ; f 表示 相机 镜头 的 焦距 ; 
民 据 式 (12) 计算 变形 后 光斑 中 心 位 置 在 物理 


D 表示 工作 距离 。 


世界 中 的 二 维 位 置 坐标 S(z。, go). 再 


通过 副 反射 面 方程 计算 变形 后 激光 光斑 的 x 坐标 ， 便 可 得 到 光斑 变形 后 的 三 维 坐标 


S(zxc, ye, ze HARAR (3) 便 可 完成 上 述 六 自由 的 计生 


5 图像 处 理 算法 实验 测试 


为 进一步 保障 文中 图 像 处 理 算 


Té E 


需 的 必要 仪器 ， 分 别 有 激 光 发 射 器 、 工 业 相 机 、 镜 头 、 移 动 平 台 以 及 激光 测 距 仪 等 。 


验 设备 的 具体 参数 如 表 1 所 示 。 


过 程 。 


YA 


行 了 相应 实验 测试 。 基 于 上 述 目 


在 本 方案 中 的 定位 精度 ， 本 文 对 


图 像 处 理 方法 的 定位 


的 以 及 第 4 章 的 图 像 整体 处 理 


流程 ， 确 定 本 实验 所 


表 1 实验 设备 参数 
设备 参数 
工业 相机 分 辨 率 : 2560x2160 (500 万 像素 ) 
HK-A5100-GM17 像 元 尺寸 : 6.5x6.5 um; 传感器 尺寸 4/3” 
工业 镜头 焦距 : 12 mm， 镜 头 尺 寸 : 1.2” 
MVL-KF1224M-25MP | 分 辩 率 : 2500 7j, C 接口 镜头 
移动 平台 移动 精度 : 10 um; 自由 度 : 单 自由 度 


实验 前 ， 首 先 将 相机 与 镜头 进行 组 装 调 试 并 预 热 ， 然后 将 激光 发 射 器 与 移动 平台 进行 组 


B. Hu 


i 


至 此 


开 激 光 器 进行 光斑 初始 位 置 的 拍摄 ， 最 后 通过 移动 平台 将 光斑 移动 


装 固定 ， 需 保障 移动 平台 与 激光 器 无 相对 运动 ， 进 而 通过 移动 平台 来 控制 光斑 移动 。 激 光 测 


距 仪 需 固 定 在 与 相机 平行 的 位 置 ， 用 以 确定 每 组 数据 的 物 距 ， 也 即 相 机 与 激光 接收 平面 的 距 
[图 6 所 示 。 KIER, 4 


E 关 闭 激 光 器 ， 利 用 相机 先 拍摄 无 光斑 时 的 背景 图 ， 然 后 打 


距离 与 真实 的 移动 距离 进行 做 差 便 可 得 到 这 一 组 实验 中 了 
物 距 以 及 移动 平台 的 移动 量 来 循环 上 述 实验 


段 已 知 的 距离 再 进行 


上 ， 一 组 完整 的 数据 便 已 采集 完成 。 按 第 4 章 所 述 对 采集 的 数据 进行 处 理 和 计算 ， 
便 可 获取 光斑 移动 前 后 的 质心 坐标 ， 进 而 得 到 了 
因为 光斑 的 位 移 是 通过 移动 平台 实现 的 ， 故 光斑 的 实际 移动 


像素 级 别 的 光斑 移动 距离 。 


E 离 是 已 知 的 。 实 验 计算 的 


两 大 组 (随机 分 ) 分 别 求 均值 和 标准 差 ， 结 果 如 表 2 所 示 。 


实验 结果 表明 ， 文 中 图 像 处 到 


E 算 法 在 处 理光 斑 图 像 时 


可 以 控制 在 较 理 想 的 亚 像素 级 别 ， 并 


的 准确 性 和 本 方案 的 稳定 性 。 
至 0.03 pixel ( 


5 像素 方法 的 定位 精度 。 通 过 改变 
操作 便 可 获得 大 量 的 结果 ， 并 把 实验 数据 分 为 


有 非常 明显 的 优势 ， 其 
两 组 数据 之 间 也 未 出 现 较 大 波动 ， 
通过 分 析 表 中 结果 ， 我 们 可 以 发 现 多 个 数据 的 定位 误差 已 经 降 
至 0.02 pixel) 的 水 平 ， 这 表明 本 文 的 图 像 处 理 算法 在 实际 应 用 中 能 够 非常 


定位 误差 
也 间接 验证 了 实验 
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6 ”实验 平台 搭建 
表 2 边缘 检测 和 质心 定位 
ibn 组 1 均值 /pixel ”组 1 标准 差 /pixel | 2H 2 均值 /pixel 组 2 标准 差 /pixel 
T y T y T y T y 

高 斯 法 0.1026 0.0347 | 0.0134 0.0191 | 0.1152 0.0606 | 0.0111 0.0208 
一 阶 重心 法 | 0.0064 0.0242 | 0.0050 0.0121 | 0.0664 0.0092 | 0.0055 0.0081 
二 阶 重心 法 | 0.0231 0.0330 | 0.0080 0.0121 | 0.0687 0.0100 | 0.0045 — 0.0059 
三 阶 重心 法 | 0.0371 0.0389 | 0.0084 0.0127 | 0.0843 0.0159 | 0.0058 0.0079 
中 心 法 0.0200 0.0586 | 0.0062 0.0148 | 0.1157 0.0387 | 0.0148 — 0.0255 
海 森 法 0.0096 0.0142 | 0.0093 0.0089 | 0.0477 0.0194 | 0.0153 0.0166 


准确 地 进行 测量 和 定位 。 另 外 ， 通 过 对 比 可 以 发 现 ， 
心 定位 。 通 过 上 述 数 据 可 以 看 到 ， 在 文中 边缘 预 处 型 
好 维持 在 0.0125 pixel 左右 。 


6 单 块 面板 的 变形 模拟 与 测量 实验 


的 加 持 下 ， 


一 阶 重心 法 更 适用 于 本 方案 的 小 目标 质 
阶 重 心算 法 的 精度 可 以 较 


为 了 验证 本 文 提出 的 面板 变形 测量 技术 的 有 效 性 和 映射 关系 的 正确 性 ， 本 文通 过 有 限 


元 方法 模拟 了 天 蕊 望远镜 主 面 的 单 块 站 


提出 的 测量 方法 , 基于 Matlab X ZRH 


板 在 重力 作用 下 的 真实 变形 ”。 


板 变形 进行 了 重 构 。 最 后 ， 通 


然后 ， 使 用 本 文 所 
过 将 有 限 元 模拟 的 真 


实 变形 结果 与 本 文 测量 结果 进行 对 比 ， 从 而 评估 该 光学 映射 测量 方法 的 误差 水 平 。 需 要 说 明 


的 是 ， 为 了 控制 单一 误差 变革 


位 误差 。 


天 马 望 远 镜 的 主 反 射 面 共 由 1008 块 实体 面板 拼接 而 成 
环 中 随机 一 块 面板 为 实例 对 象 ， 并 所 
后 模拟 重力 变形 后 的 坐标 ， 作 为 真实 上 
行 实验 ， 并 与 真实 变形 进行 对 比 。 首 # 


二 和 科学 验证 六 自 1 


度 映 射 测量 法 ， 因 


， 共 分 


此 光斑 位 移 量 并 未 引入 定 
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E 标 是 已 知 的 。 


区 了 该 面板 在 俯仰 角 分 别 为 45° 和 15° 时 的 有 限 元 ， 然 
的 面板 变形 情况 。 随 后 ， 采 / 
Et， 面板 变形 前 4 角 点 的 4 


j 本 文 所 提出 的 测量 方法 进 


通过 面板 角 点 
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标 与 天 线 焦点 坐标 可 以 得 到 各 激光 光束 的 直线 方程 ， 并 联 立 天 线 副 面 方程 求解 激光 与 副 
面 的 交点 坐标 ， 也 即 光 斑 坐 标 。 用 同样 方法 ， 可 以 求 得 变形 后 光斑 的 位 置 坐标 ， 从 而 获得 准 
确 的 光斑 移动 量 。 接 下 来 ， 利 用 光斑 的 移动 量 重 构 面 板 的 变形 ， 即 第 3 章 求解 刚体 面板 的 
空间 六 自由 度 。 已 知 面板 各 角 点 的 六 自由 度 变 形 量 后 ， 基 于 变形 前 各 角 点 的 原始 坐标 ， 可 以 
得 到 面板 变形 后 的 各 角 点 坐标 ， 这 些 坐 标 即 为 本 文 方法 的 测量 结果 。 最 后 ， 将 测量 结果 与 有 


限 元 模拟 变形 后 的 坐标 进行 比较 ， 得 出 本 文 方法 的 测量 误差 。 详 细 结果 数据 见 表 3 和 表 4。 


[之 


» 


id) 


S 


表 3 俯仰 角 45 时 重力 变形 的 测试 结果 
点 编号 
参量 519 520 615 616 
X | —16.58084588 —17.37749485  —18.23593996 —19.112 10009 
Ru Y | 12.7215286 11.60985744 13.991 41831 12.768 777 41 
Xm Z | 5.51740581 5.251690669 6.352 309 82 6.352 270 594 
X 0.39 0.391 0.503 0.505 
光斑 位 移 量 /mm | Y 1.71 1.787 1.629 1.718 
Z 0.304 0.277 0.269 0.243 
X | —16.5808673  —17.37748224  —18.23594944 — —19.1120818 
本 文 测量 结果 /m | Y | 12.72152635 11.609 835 61 13.991 439 28 12.768 780 52 
Y | 5.251743866 5.251683756 6.352325076 6.352258966 
X | 2414x107  —126x10 9? 9.48x 107° —1.83 x 107% 
误差 /m Y | 225x10 96 2.18 x 107°%  —240x10 9 一 3.11 x 10-96 
Z | -—329x10 99  Á 691x106  Á —153x10 9 1.16 x 107% 
FEI /s 1.364 657 981 
表 4 俯仰 角 15? 时 重力 变形 的 测试 结果 
点 编号 
参量 519 520 615 616 
FONS X | —16.58040933  —17.37721667  —18.23546282  —19.11171886 
dud Y | 12.715087 09 11.60351509  13.98446446 12.761906 58 
Z | 562280514 5.261890323 . 60.363400155 6.363 126 344 
X 4.388 4.468 5.191 5.27 
光斑 位 移 量 /mm | Y 9.202 9.787 8.62 9.286 
Z 0.87 0.703 0.502 0.341 
X | —16.580 420993  —17.377142 76 |—18.23550228  —19.111636 69 
本 文 测量 结果 /m | Y | 12.71508125 11.603502772 13.984475877 12.76192142 
Y | 5.262278444 5.261870998 . 6.363498005 6.363 049 888 
X 1.17 x 10795 —7.39x10 9 3.95 x 107% 8.22 x 107% 
误差 /m Y 5.84 x 10-96 1.23 x 107%  —114x10 9  —148x10 95 
Z 2.07 x 10-96 1.93 x 107° . —9.79x10 9 — 7.65 x 107% 
耗 时 /s 1.548967 489 
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X 3 和 表 4 中 提供 了 有 限 元 模拟 重力 变形 后 的 真实 角 点 4 
动量 以 及 变形 测量 结果 、 测 量 误差 、 耗 
刚体 面板 姿态 变化 的 测量 结果 达到 了 预期 的 效果 。 与 有 限 元 仿 
微米 到 几 十 微米 的 范围 内 。 这 一 结果 进一步 证 实 了 本 文 所 提出 的 光学 映射 测量 法 能 够 准确 
地 测量 大 型 射电 望远镜 天 线 面板 在 六 个 自由 度 上 的 变形 ， 包 括 平 移 、 旋 转 等 多 种 变形 形式 。 
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标 和 本 文 方法 得 到 的 光斑 移 
时 数据 。 从 误差 结果 可 知 ， 利 用 本 文 映 射 方法 对 天 线 
真 结果 相 比 ， 测 量 误差 仅 在 几 


此 外 ， 利 用 本 文 方法 ， 安 装 好 测量 系统 后 ， 单 块 面板 的 变形 仅 在 2 s 内 便 可 快速 完成 重 构 ， 


这 对 于 天 线 局 部 测量 与 评估 具有 重要 意义 ， 并 为 实时 监督 和 调整 天 线 面板 的 姿态 提供 了 可 


AERE BR. 


7 fh i8 


本 文 提出 了 
度 姿态 变形 理 


全 新 的 基于 光学 方法 的 天 线 面 板 变形 习 
论 对 面板 刚性 变形 进行 表征 ， 将 激光 光斑 与 天 线 面 板 变 形 进行 关联 ， 并 成 功 


E 构 技术 。 该 技术 采用 刚体 六 自 


E 导 出 了 两 者 之 间 的 映射 关系 。 将 这 一 映射 关系 与 六 自 1 


Z] 


行 了 改进 。 本 文 将 


光斑 点 这 样 的 小 目标 时 容易 受到 背景 干扰 而 导致 识别 准确 


| 


Hr 


y» 


别 的 定位 ， 进 


未 来 还 可 以 进一步 优化 本 文 方法 ， 例 如 解 耦 天 线 所 有 面板 之 间 


至 天 线 所 有 面板 上 ， 对 天 线 面 型 进行 实时 的 测量 。 另 外 


和 速度 ， 并 看 


完 各 种 误差 的 补偿 方法 等 ， 扩 大 应 用 范围 


度 变 形 


性 下 
像 差 分 方法 与 Canny 边缘 检测 算法 相 结合 ， 
别 的 同时 ， 有 效 去 除 背 景 干扰 ， 提 高 光斑 质心 的 准确 性 和 稳定 性 。 另 外 ， 为 了 进一步 提高 》 
斑 质 心 的 定位 精度 ， 本 文 还 采用 了 高 精度 亚 像素 质心 定位 
提高 了 整个 光斑 图 像 处 理 算法 的 准确 性 和 效率 。 


I 


Ux 


景 ， 为 天 线 领域 的 发 展 做 出 更 多 的 贡献 。 
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Research on Deformation Reconstruction Technology of 


Antenna Panel Based on Laser Mapping Method 


YANG De-jin!l2, ZHANG Zihan?, YE Qian?4 
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Tong University, School of Mechanical and Power Engineering, Shanghai 200240, China; 4. Key 
Laboratory of Radio Astronomy, Chinese Academy of Sciences, Nanjing 210008, China) 


Abstract: In this paper, a single panel deformation of the main reflector of a large radio 
telescope is studied and the six degrees of freedom attitude change of robotics is applied to 
the deformation of the reflector panel of a large radio telescope. A characterization method 
of panel deformation based on six degrees of freedom deformation of rigid body is proposed. 
Based on the six-DOF attitude deformation theory of rigid body, an optical reconstruction 
technique of antenna panel deformation is proposed in this paper. This technology can 
reconstruct the deformation of the panel quickly and accurately by the spot image presented 
by the laser unit consolidated on the panel. In this paper, the mapping relationship between 
laser spot and antenna panel deformation is deeply explored, and the spatial deformation 
problem of antenna panel is transformed into the laser spot displacement problem in a finite 
region. And the light spot image of this scheme is studied, and an effective light spot image 
processing algorithm is provided. It integrates the image difference method and Canny edge 
detection algorithm, which can effectively solve the problem of traditional Canny algorithm 
in dealing with small targets. It also uses the high-precision sub-pixel centroid localization 
algorithm, which can locate light spots at sub-pixel level. It provides a new and reliable 
image processing method for antenna panel measurement. This study also provides a new 


idea and method for improving the accuracy and efficiency of antenna panel measurement. 


Key words: large radio telescopes; deformation reconstruction; image processing; subpixel 


